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Abstract 

Genetic differentiation among population samples of G. 
gauthieri on both sides of the Strait of Gibraltar was investigated 
by enzyme electrophoresis at 20 enzyme loci with reference to G. 
gauthieri and G. ibericus from different areas in the Iberian 
peninsula. 

Levels of divergence resolved suggest that G. gauthieri 
reached North Africa from the Iberian peninsula, and not vice 
versa. With reference to geological evidence, this was likely to 
occur during the Messinian (6-5 My ago). As the opening of the 
Strait of Gibraltar represents a well dated geological event, an 
attempt has been made to calibrate the molecular clock for the 
G. pule: c-group. 

Indirect estimations of gene flow levels were shown to be ex¬ 
tremely low. Probably as a consequence of frequently occurring 
semi-arid conditions, populations of G. gauthieri are strongly 
subdivided and likely to be subject to stochastic processes. 
Therefore, even within areas, population samples may be geneti¬ 
cally markedly different. 

Dendrograms of genetic distances were compared with geo¬ 
logical evidence. From these data an hypothesis of phylogenetic 
relationships is proposed. 

Resume 

La differentiation genetique entre des pr£lfevements de popula¬ 
tions de G. gauthieri de chaque cote du detroit de Gibraltar, par 
rapport a G. gauthieri et G. ibericus issus de diff£rentes regions 
de la Peninsule lb£rique, a et£ Studiee par £lectrophor£se a 20 
loci d’enzymes. 

Les niveaux de divergence observes suggerent que G. gauthieri 
a gagn£ l’Afrique du Nord a partir de la Peninsule Ib6rique, et 
non pas vice versa. Compte tenu des donnees geologiques, cela 
se serait produit pendant le Messinien (il y a 6-5 Ma). Puisque 
l’ouverture du detroit de Gibraltar represente un 6venement g6- 
ologique bien date, un essai a ete fait de calibrer l’horloge 


moleculaire pour le groupe G. pulex. 

Des estimations indirectes des niveaux de flux genetique se 
sont montrees extremement basses. Probablement a la suite des 
conditions de semi-aridite intervenant frequemment, les popula¬ 
tions de G. gauthieri sont fortement subdivisees et susceptibles 
d’etre soumises k des processus stochastiques. De ce fait, des 
prelevements de populations peuvent etre genetiqement remar- 
quablement differents, meme a Tinterieur d’une meme region. 

Des dendrogrammes de distances genetiques ont etes mises en 
rapport avec des faits geologiques. A partir de ces donnees, une 
hypothese des relations phylogenetiques est proposee. 


Introduction 

The Iberian members of the Gammarus pulex - 
group (sensu Karaman & Pinkster, 1977) G . 
gauthieri S. Karaman, 1935 and G. ibericus Mar- 
galef, 1951 were shown to be genetically very dis¬ 
tinct from taxa of this group from north of the 
Pyrenees. According to Scheepmaker et al., (1988), 
with reference to other species, genetic differentia¬ 
tion between French and Iberian species was ob¬ 
served at the generic level. As the freshwater fauna 
of the Iberian peninsula is related to that of North 
Africa, Scheepmaker et al. (1988) suggested the 
possibility that Iberian members of the G . pulex- 
group originate from forms in North Africa, and 
not from the closely related forms from western Eu¬ 
rope north of the Pyrenees. Although G. gauthieri 
is distributed in Spain, Morocco, Algeria and 
Tunisia, it is not clear whether it reached Morocco 
from the Iberian peninsula or vice versa (Scheep¬ 
maker et al., 1988). G. ibericus is endemic for the 
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Table L Sampling localities and species distribution. 


Station 

no. 

Area in 
fig. 1A, B 

Species 

Country 

Prov /dept 

Drainage 

system 

Locality description/reference 

Sampling 

date 

1 

- 

G. d. duebeni 

Holland 

North-Holland 

• 

Brackish pool behind the dikes of the Hondsbosche 
Zeeweering at Petten, westside of the Vereenigde 

Harger- and Pettemerpolder, near road N9 

3-VII-*88 

2 

- 

G . d. duebeni 

Fiance 

Pas-de-Calais 

* 

Rockpool along the coast at Audresselles, near 
road D940 

24-VU-'88 

3 

C 

G. d. celticus 

France 

Finistfcre 

LaFlfcche 

At the crossing with road D788,1.5 km E. of Lesneven 

25-Vn-"88 

4 

C 

G . d. celticus 

France 

Finistfcre 

Elom 

Tributary, at the municipal camping of Sizun, near 
road D764 

25-VU-'88 

5 

* 

G. gauthieri 

Morocco 

Marrakech 

OuedAouri 

Source at Imi-n-Tala, at the end of a mountain track 

7 km S. of Amizmiz and road S507 

23-IIK88 

6 

D 

G. gauthieri 

Morocco 

Khenifra 

- 

Stream draining directly in the Aguelmane de Sidi AH, 
near road 3211,20 km S. of Timhadite 

23-IIK88 

7 

D 

G. gauthieri 

Morocco 

Khenifra 

Oued Beth 

Ain Achbal, source of the Tigrigra tributary, 2 km W. of 
Azrou, near road P24 

24-IIK88 

8 

D 

G. gauthieri 

Morocco 

Khenifra 

Oued Sebou 

Oued Mikkfcs, tributary alongside road S309,3 km S.E. 
of Ifrane 

24-III- , 88 

9 

D 

G. gauthieri 

Morocco 

Khenifra 

Oued Sebou 

Tributary at the crossing with road 4629,10 km S.S.E 
of Imouzzfcr-du-Kandar 

24-III-88 

10 

E 

G. gauthieri 

Spain 

Ronda 

Rio Guadiaro 

Rio Guadalobacin, tributary crossing road 420, 

1.5 km N. of Ariate 

26-IIK88 

11 

E 

G. gauthieri 

Spain 

Ronda 

Rio Guadalhorce 

Rio de las Cuevas, near road C341,0.5 km E. of 

Cuevas del Becerra 

27-m-'88 

a 

F 

G. gauthieri 

Spain 

Cuenca 

Tajo 

Station 6 in Scheepmaker et al. (1988) 

- 

b 

F 

G. gauthieri 

Spain 

Cuenca 

Guadiana 

Station 7 in Scheepmaker et al. (1988) 


c 

G 

G. gauthieri 

Spain 

Valladolid 

Duero 

Station 10 in Scheepmaker et al. (1988) 

- 

d 

G 

G. gauthieri 

Spain 

Valladolid 

Duero 

Station 11 in Scheepmaker et al. (1988) 


e 

G 

G. gauthieri 

Spain 

Valladolid 

Duero 

Station 12 in Scheepmaker et al. (1988) 

- 

f 

F 

G. ibericus 

Spain 

Cuenca 

Jucar 

Station 8 in Scheepmaker et al. (1988) 


g 

F 

G. ibericus 

Spain 

Cuenca 

Tajo 

Station 9 in Scheepmaker et al. (1988) 



Serrama de Cuenca, a mountainous area in central 
Spain (Lop, 1987). 

The present study investigates the possible ori¬ 
gins of G. gauthieri and G. ibericus. Primarily, 
genetic relationships among G. gauthieri popula¬ 
tions from North Africa and the Iberian peninsula 
on either side of the Strait of Gibraltar are studied 
with reference to levels of divergence among popu¬ 
lations of G. gauthieri and G. ibericus resolved 
within the Iberian peninsula by Scheepmaker et al. 
(1988). Intra-specific genetic differentiation was ex¬ 
amined in relation to estimated levels of gene flow. 
As indices of gene flow based on F-statistics may be 
influenced by natural selection and therefore may 
constitute poor indicators of the genetic structure 
of populations, these indices were used in combina¬ 
tion with the methods of Slatkin (1981, 1985); the 
latter methods are believed to be less sensitive to the 
influence of natural selection (Barton & Slatkin, 
1987). 

On the other hand, divergence data and phylo¬ 


genetic relationships are compared with geological 
evidence. The most obvious and recent migratory 
pathway between the Iberian peninsula and North 
Africa is across the Strait of Gibraltar (e.g. Busack, 
1986, Scheepmaker et al., 1988). Extension of G. 
gauthieri may have occurred during the Messinian 
salinity crisis, when the Mediterranean dried up 
(6-5 My ago; e.g. Hsu et al., 1977). However, the 
complexity of the paleogeography of the western 
Mediterranean (e.g. Dercourt et al., 1985; Mal¬ 
donado, 1985) is an obstacle to any a priori assump¬ 
tion concerning the distributional history of G. 
gauthieri. Therefore, other hypotheses of dispersal 
or vicariant events leading to its present distribu¬ 
tion should not be discarded. 

It may be justified to calibrate the molecular 
clock for a given set of loci in a group of closely 
related species, as long as the evolutionary rates for 
the individual loci remain constant (Scheepmaker et 
al., 1988). As the results of this study are in favour 
of a Messinian extension across the Strait of Gibral- 
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Fig. 1A-E . Study areas and sampling localities (cf. Table I). 


tar, an attempt was made to provide an initial 
calibration of the molecular clock for the G. 
pu lex-group. 

Genetic variation at 20 presumptive enzyme loci 
was investigated using starch gel electrophore¬ 
sis. A network of phylogenetic relationships was 
rooted by the outgroup method (Farris, 1972). 
In addition to Scheepmaker et al. (1988), hybridi¬ 
zation experiments among samples of G. gauthieri 
from Morocco and the Iberian peninsula were 
carried out to evaluate the taxonomic status of 
the forms involved, and the evolutionary sig¬ 
nificance of the genetic variation recorded. From 
these data, an evolution model of G. gauthieri 
in the western Mediterranean area is proposed. 


Materials and methods 

Sampling and collection sites 

Sampling for electrophoretic studies and cross¬ 
breeding experiments was carried out based on the 
procedures of Scheepmaker (1987). The collection 
sites and sampling dates along with the species col¬ 
lected are listed in Table I and Figs. 1A-E. In 
Morocco, the area closest to the Strait of Gibraltar 
from which gammarids were reported, was sampled 
at four localities to evaluate local genetic variation 
(the Khenifra area: Fig. 1A, area D; Table I and 
Fig. ID, stations 6-9). One of these localities (sta¬ 
tion 6) is the type locality of G. acalceolatus Pink¬ 
ster, 1971 (a form morphologically very similar to 
G. ibericus , but calceoli are lacking). However, 
only G. gauthieri was sampled at this locality. To 
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evaluate geographic variation over greater dis¬ 
tances, one sample in the High Atlas was taken 
(Table I, Fig. 1A; station 5). In Spain, the area 
closest to the Strait of Gibraltar from which gam- 
marids are known is the province Ronda (Table I, 
Figs. 1A, E, area E; stations 10-11). In that area, 
due to an increased pollution of most tributaries, 
not more than two localities could be sampled. For 
purposes of outgroup comparison, two samples of 
the marine and brackish form G . duebeni duebeni 
Liljeborg, 1852 (Table I, Fig. 1A; stations 1-2), 
and two samples of the freshwater form G. duebeni 
celticus Stock & Pinkster, 1970 (Table I, Fig. IB, C, 
area C; stations 3-4) were studied simultaneously. 


Electrophoresis 

Electrophoresis and staining procedures were iden¬ 
tical to those of Scheepmaker et al. (1988). The fol¬ 
lowing enzyme systems were assayed (E. C. num¬ 
bers between parentheses): ADA - Adenosine 
deaminase (3.5.4.4); ALP-Alkaline phosphatase 
(3.1.3.1); APK-Arginine phosphate kinase 
(2.7.3.3); EST-Esterase (3.1.1.1); GDH-Glu¬ 
tamate dehydrogenase (1.1.1.47); GOT - Glu¬ 
tamic oxaloacetic transaminase (2.6.1.1); GPI- 
Glucose phosphate isomerase (5.3.1.11); HK- 
Hexokinase (2.7.1.1); LAP - Leucine aminopep- 
tidase (3.4.1.1); MDH - Malate dehydrogenase 
(1.1.1.37); ME-Malic enzyme (1.1.1.40); MPI- 
Mannose phosphate isomerase (5.3.1.8); PEP- 
Peptidase (3.4.11/13); 6PGD-6 Phosphoglu- 
conate dehydrogenase (1.1.1.44); PK-Pyruvate 
kinase (2.7.1.40). 


Analysis of allozyme variation 

The genetic interpretation of the variation was in¬ 
ferential. Data were analyzed by the computer pro¬ 
gram BIOSYS-1 (Swofford & Selander, 1981). A 
matrix of Nei’s (1972) genetic distance (D) was cal¬ 
culated for all populations sampled for this study, 
and for G. gauthieri and G. ibericus from the Iberi¬ 
an peninsula sampled by Scheepmaker et al. (1988). 
From these data (for sake of space not all presented 


here), an UPGMA dendrogram (Sneath & Sokal, 
1973) was constructed. A “character state” trans¬ 
formation scheme was generated with the aid of the 
program Jelly (Ellis, 1987). The Jelly network 
generates ancestors (Hypothetical Taxonomic 
Units or HTU’s) with their own gene frequencies 
(thus not selected from one of the states occurring 
in the Operational Taxonomic Units or OTU’s) that 
sum up to one using Rogers’ (1984) HAP al¬ 
gorithm; for a discussion of this method see 
Scheepmaker et al. (1988). 

The UPGMA method and Nei’s (1972) distance 
coefficient require both the assumption of homo¬ 
geneous evolutionary rates (e.g. Farris, 1981). 
Moreover, Nei’s distance is unsuitable for deriving 
phylogenies because of its non-metricity. The Jelly 
network does not require any assumption about 
evolutionary rates. Rooted by the outgroup 
method, this procedure seems more appropriate to 
infer phylogenetic relationships; however, the 
choice of an outgroup taxon is not directly obvious. 
Although genetically distant and reproductively 
isolated from G. gauthieri and G. ibericus (Scheep¬ 
maker et al., 1988), members of the G. pulex- group 
in northwestern Europe are genetically less dif¬ 
ferentiated, suggesting a more recent origin. There¬ 
fore, they may not represent a suitable outgroup. 
On the other hand, related forms may be genetically 
so dissimilar, that they do not contribute to resolve 
phylogenetic relationships (e.g. Skadsheim & Sie- 
gismund, 1986). Genetically as well as morphologi¬ 
cally, G. duebeni (i.e. G. d. duebeni + G. d . celti¬ 
cus) appears to be sufficiently distinct from the in¬ 
group to meet the requirements of an outgroup 
without being too different. The state of several 
taxonomically important characters in G. duebeni 
is intermediate between that encountered in the G. 
pulex-g roup and in the systematically distinct 
brackish and marine forms like G. zaddachi Sex¬ 
ton, 1912 and G. salinus Spooner, 1947 (viz. pro¬ 
truding lobe backwards of basipodite P5-7, 
setation and armature P3-7, setal fringe of the 
mandibular palp, and shape of the eye; Karaman & 
Pinkster, 1977; Lincoln, 1979; Gledhill et al., 1976; 
Pinkster & Platvoet, 1986). G. duebeni is common¬ 
ly assumed to be old and formerly widespread in 
western Europe (e.g. Hynes, 1954, 1955; Pinkster 
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Table II. Electromorph frequencies at 17 presumptive gene loci. 
d ~ G. d. duebeni ; c - G. d. celticus ; g = G. gauthieri; 
N = sample size; h = heterozygosity per locus (direct count); 
H - mean heterozygosity over all (including monomorphic) 
loci; * = significant departure from Hardy-Weinberg distribu¬ 
tion (P < 0.05); electromorph designation from fastest to 
slowest in alphabetical order, in accordance with Scheepmaker 
& Van Dalfsen (1989) and Scheepmaker (1989); alternative 
allelic designations employed by Scheepmaker et al. (1988) are 
given in parentheses. 


staiion 1 2 3 4 5 6 7 8 9 10 11 

ddccggggggg 

locus 

electromorph 


Alp-l 

(N) 

52 

8 

13 

20 

20 

20 

20 

30 

30 

30 

30 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 


b 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

Alp-2 

(N) 

52 

15 

20 

30 

20 

20 

20 

30 

30 

30 

30 

a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

(a) 

c 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

Apk 

(N) 

63 

25 

30 

20 

25 

30 

29 

40 

35* 

25 

35 

a 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

0.40 

1.00 

1.00 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.60 

0.00 

0.00 


c 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 




- 

' 


- 

• 

0. 23 



Est-2 

(N) 

55 

25 

30 

35 

30 

40 

20 

25 

33 

35 

35 


a 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(a) 

b 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

Gdh 

(N) 

53 

25 

30 

35 

40 

40 

30 

45 

48 

43 

45 


a 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

0.00 

0.00 

1.00 

1.00 


b 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

Got-1 

(N) 

55 

31 

36 

25 

46 

43 

46 

41 

51 

41 

35 

(a) 

b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.89 

0.00 

0.00 

(b) 

c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.06 

0.97 

0.01 

0.00 

0.94 

0.80 

e 

0.00 

0.15 

0.00 

0.36 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.04 


f 

0.00 

0.00 

0.00 

0.00 

1.00 

0.94 

0.03 

0.92 

0.00 

0.06 

0.16 


g 

0.11 

0.11 

1.00 

0.64 

0.00 

0.00 

0.00 

0.07 

0.00 

0.00 

0.00 


h 

0.88 

0.74 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.11 

0.00 

0.00 


i 

0.01 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.16 

0.36 


0.40 

- 

0. 07 

0. 07 

0.17 

0. 21 

0.12 

0.31 

Got-2 

(N) 

54 

31 

36 

25 

46 

43 

46 

40 

51 

41 

35 

(a) 

b 

0.01 

0.00 

0.00 

0.00 

0.00 

0.98 

1.00 

1.00 

1.00 

1.00 

1.00 

c 

0.99 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(O 

d 

0.00 

0.00 

0.00 

0.00 

1.00 

0.02 

0.00 

0.00 

0.00 

0.00 

0.00 

h 

0.02 




- 

0.05 

- 

• 

- 

* 

- 

Gpi 

(N) 

87 

39 

47 

53* 

55 

43 

27 

52 

53 

45 

39 

e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.07 

0.00 

0.04 

0.00 

0.00 

0.00 


f 

0.00 

0.00 

0.00 

0.00 

0.00 

0.85 

0.04 

0.00 

0.00 

0.46 

0.12 


g 

0.00 

0.00 

0.00 

0.00 

0.96 

0.01 

0.70 

0.81 

0.76 

0.00 

0.00 


h 

0.00 

0.00 

0.00 

0,00 

0.00 

0.00 

0.06 

0.00 

0.02 

0.39 

0.56 


i 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.20 

0.00 

0.00 

0.03 

0.05 

(h) 

m 

0.01 

0.03 

0.00 

0.00 

0.00 

0.00 

0.00 

0.15 

0.22 

0.12 

0.27 

0) 

n 

0.98 

0.91 

0.15 

0.48 

0.00 

0.07 

0.00 

0.00 

0.00 

0.00 

0.00 

0) 

0 

0.01 

0.06 

0.85 

0.52 

0.05 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

h 

0.03 

0.18 

0.21 

0. 28 

0.09 

0.26 

0. 59 

0.03 

0.44 

0.53 

0.51 

Hk 

(N) 

41 

26 

26 

30 

61 

33 

36 

46 

39 

51 

45 


b 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


c 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Mdh-1 

(N) 

49 

41 

36 

25 

46 

40 

46 

41 

51 

41 

31 


e 

0.00 

0.00 

0.00 

0.04 

0.00 

0.00 

0.00 

0.00 

0.34 

0.00 

0.00 

(c) 

f 

0.94 

0.99 

0.00 

0.08 

1.00 

0.00 

0.94 

0.99 

0.04 

1.00 

1.00 

(d) 

g 

0.06 

0.01 

1.00 

0.88 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(e) 

h 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.07 

0.01 

0.62 

0.00 

0.00 

i 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.08 

0.02 


0. 24 


- 

0. 09 

0. 02 

0.49 

• 


Mdh-2 

(N) 

47 

41 

36 

25 

38 

43 

46 

41 

51 

41 

35 

(a) 

b 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

(b) 

c 

0.00 

0.00 

0.00 

0.00 

0.32 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

d 

0.00 

0.00 

0.00 

0.00 

0.68 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 




- 

0.32 

- 

- 

- 


- 

* 


35 


Me 

(N) 

52 

31 

36 

38 

55 

52 

44 

54 

62 

49 

51 


a 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.00 

1.00 

1.00 


b 

0.01 

0.00 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 


c 

0.99 

1.00 

0.00 

0.00 

0.00 

0,00 

0.00 

0.00 

0.00 

0.00 

0.00 


d 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.00 

0.00 


h 

0.02 







- 



* 

Mpi 

(N) 

66 

31 

28 

33 

3 

23 

31 

27 

35* 

31 

30 

f 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(O 

g 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

h 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

0.98 

0.83 

0.01 

0.00 

0.00 


i 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.17 

0.90 

0.00 

0.00 

(0 

k 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.09 

0.00 

0.00 

h 





- 


0.03 

0.11 

0.16 


- 

Pep-2 

(N) 

51 

34 

31 

33 

31 

38 

36 

22 

37 

31 

20 

c 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.02 

0.00 

1.00 

1.00 


d 

1.00 

1.00 

1.00 

1.00 

0.98 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(d> 

e 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

0.98 

0.88 

0.00 

0.00 

g 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.00 

0.00 

0.00 

0.00 

0.02 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 


j 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.12 

0.00 

0.00 


h 

- 




0.03 

* 


0.04 

0. 24 



Pep-4 

(N) 

31 

34 

31 

33 

31 

36 

36 

21 

37 

31 

20 


b 

0.00 

0.00 

0.00 

0.00 

0.00 

0.92 

1.00 

1.00 

0.89 

1.00 

1.00 

(a) 

c 

0.48 

0.47 

0.00 

0.38 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

(c) 

d 

0.32 

0.12 

0.00 

0.00 

1.00 

0.08 

0.00 

0.00 

0.11 

0.00 

0.00 


g 

0.19 

0.24 

1.00 

0.62 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.00 

0.18 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


h 

0.05 

0.64 


0.46 

- 

0.11 

- 


0.22 



6-Pgd 

(N) 

30 

25 

30 

30 

20 

10 

20 

20 

20 

20 

20 


b 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


c 

1.00 

1.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Pk 

(N) 

56 

51 

51 

48 

36 

51 

56 

41 

47 

51 

40 


a 

0.00 

0.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 


c 

0.00 

0.00 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


e 

h 

1.00 

1.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


H 

0.04 

0.06 

0.01 

0. 07 

0.02 

0.02 

0.04 

0.03 

0.10 

0.03 

0.04 


et al., 1970; Golikov & Tzvetkova, 1972); com¬ 
pared with members of the G. pulex-g roup as well 
as with G. zaddachi and G. salinus, reproduction 
occurs at a slow rate which is considered a primitive 
character by Pinkster (pers. comm.). 

To measure population subdivision, F-statistics 
(0) were calculated according to Weir & Cockerham 
(1984). For loci with multiple alleles, a weighted 
average 8 was calculated. Significance of the ob¬ 
served 0-values was tested with a chi-square heter¬ 
ogeneity test (Workman & Niswander, 1970). Gene 
flow estimates (Nm) were derived from the relation¬ 
ship^ 1/1 +4 Nm, where TVis the effective popula¬ 
tion size and m the migration rate (Wright, 1931). 
8- and Ato^-values were calculated with the aid of 
the program Theta (Ellis, 1989). 

Gene flow levels (Nm) among populations and 
areas were also calculated from the average fre¬ 
quency of private alleles p (1) according to Slatkin 
(1985). To adjust for differences in sample size, a 
corrected Nm-value, Nm*, is obtained by dividing 
the Nm estimate by the ratio of the actual sample 
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Table III. Number of diagnostic loci (0.99 criterion; Ayala & Powell, 1972) among samples of the populations studied, averaged by 
area. Data from Scheepmaker et ah (1988) are included. 



■re* / 
llalian(s) 

no. of 
samples 

1.2 

C; 3,4 

5 

D;6 

D-.7.8.9 

E; 10.11 

F;a.b 

G; c,d,c 

F; f 

F;g 

G. <L duebeni 

1.2 

2 

- 










G. d. cebicus 

C;3,4 

2 

2 

- 









G. gauthieri 

5 

1 

10 

12 

- 








G. gauthieri 

D;6 

1 

11 

14 

• 

- 







G. gauthieri 

D-,7,8.9 

3 

9 

13 

3 

3 

• 






G. gauthieri 

E; 10.11 

2 

11 

13 

8 

7 

4 

- 





G. gauthieri 

F;«.b 

2 

12 

14 

7 

5 

3 

5 

• 




G. gauthieri 

G; c,d,e 

3 

11 

12 

9 

6 

2 

8 

9 

- 



G. ibericus 

F; f 

> 

12 

12 

6 

5 

2 

6 

3 

6 

- 


G. ibericus 

F;g 

1 

10 

11 

10 

10 

7 

11 

8 

9 

5 

- 


size to 25 (Slatkin, 1985). In addition, gene flow 
profiles were obtained by plotting the values of the 
conditional average frequencies (CAF) according 
to Slatkin (1981). These procedures have been ap¬ 
plied to the complete data set and to different sub¬ 
sets of populations. For a more detailed discussion 
of these methods and their application within the 
G. pulex- group is referred to Scheepmaker (1990). 

Cross-breeding experiments 

Cross-breeding experiments were conducted ac¬ 
cording to Pinkster (1983) and were carried out (a) 
between geographically close population samples 
of G. gauthieri from Morocco (station 6, area D) 
and southern Spain (station 10, area E); (b) between 
population samples of G . gauthieri from Morocco 
(station 6) and central Spain (station b, area F) and 
(c) between population samples of G. gauthieri 
from Morocco (station 6) and G. ibericus (station 
g, area F). 


Results 

Twenty presumptive gene loci from 17 enzyme sys¬ 
tems were scorable in all populations. Electro- 
morph frequencies and relative mobilities are listed 
in Table II. The Ada , Est-l and Lap loci were 
monomorphic and fixed for the same electromorph 
in all population samples. The number of loci and 
subunit structure were consistent with the interpre¬ 
tation of Scheepmaker et al. (1988). 

Genetic variation among populations 

Differences in electromorph distribution among the 
samples studied are shown in Table II. G. gauthieri 
population samples from areas D, E and station 5 
are characterized by the Apk °, Est-2 b , 6Pgd b , Pk° 
and, except for station 5, by the Got-2 b and Pep-4 b 
electromorphs. 

The Alp-l a , Alp-2 a > Mpi g and Pep-2 C electro¬ 
morphs are restricted to G. gauthieri samples from 
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Table IV. Matrix of Nei’s (1972) genetic distance coefficients among samples of the populations studied, averaged by area. Data from 
Scheepmaker et al. (1988) are included. Range D mjn - D max in italics; (*) = only one population sample included. 


species 

ares/ 

stilion(s) 

no. of 
samples 

1.2 

C:3.4 

5 

D; 6 

D; 7,8,9 

E; 10,11 

F; a,b 

0;c,d,e 

F; f 

G . d . duebeni 

1.2 

2 

0.00 

(0.00-0.00) 









G . d . celticus 

C; 3,4 

2 

0.25 

(0.21-0.29) 

0.02 

(0.02-0.02) 








G . gauthieri 

5 

1 

0.98 

(0.98-0.99) 

1.16 

(1.13-1.18) 

C) 







G . gauthieri 

D; 6 

1 

1.32 

(1.32-1.32) 

1.34 

(1.32-1363) 

0.30 

(030-030) 

C) 






G . gauthieri 

D-.7.8.9 

3 

1.09 

(1.00-1.18) 

1.20 

(1.10-1.36) 

0.41 

(033-036) 

032 

(0.25-0.42) 

0.20 

(0.10-030) 





G . gauthieri 

E; 10.11 

2 

1.34 

(1.33-1.35) 

135 

(132-1.58) 

0.65 . 
(0.64-0.65) 

0.54 

(033-0.54) 

0.49 

(033-0.68) 

0.01 

(0.01-0.01) 




G . gauthieri 

F; a,b 

2 

1.62 

(1.53-1.70) 

1.96 

(1.84-2.04) 

0.61 

(035-0.68) 

0.47 

(0.41-032) 

0.47 

(033-0.59) 

031 

(030-033) 

0.10 

(0.10-0.10) 



G . gauthieri 

G; c.d.e 

3 

0.95 

(0.92-0.98) 

1.08 

(1.02-1.14) 

0.79 

(0.71-0.85) 

0.62 

(036-0.66) 

0.47 

(039-0.51) 

0.74 

(0.67-0.79) 

0.70 

(0.60-0.78) 

0.04 

(0.02-0.05) 


G . ibericus 

F; f 

1 

1.10 

(1.01-1.11) 

1.32 

(1.26-1.39) 

0.48 

(0.48-0.48) 

0.48 

(0.48-0.48) 

038 

(0.290.43) 

0.47 

(0.47-0.47) 

0.23 

(0.22-0.24) 

0.55 

(030-0.57) 

<*> 

G . ibericus 

F;g 

1 

1.22 

(1.21-1.23) 

1.48 

(1.41-136) 

0.79 

(0.79-0.79) 

0.82 

(0.82-0.82) 

0.74 

(0.71-0.79) 

0.93 

(0.91-0.94) 

0.66 

(0.65-0.66) 

1.06 

(0.98-1.09) 

039 

(039-039) 


area E; the Pep-2 d electromorph characterizes G. 
gauthieri samples from stations 7, 8 and 9 in area D. 

Electromorphs of other polymorphic loci are ir¬ 
regularly distributed within and among areas. 
Within area D, at the Me locus, G. gauthieri sam¬ 
ples of stations 6, 7 + 8, and 9 are fixed for three 
different electromorphs: Me b , Me a and Me e , 
respectively; Me a occurs also in G. gauthieri sam¬ 
ples from area E, stations 10-11; Me b is also 
found in the G. gauthieri sample from station 5. 
Likewise, in area D, G. gauthieri samples from sta¬ 
tions 6 and 9 are characterized by the Mdh-V and 
Mdh-l e > h electromorphs, that do not occur in 
other samples; at the Gdh locus, G. gauthieri sam¬ 
ples from stations 8 and 9 are, unlike station 7, 
fixed for the Gdh b electromorph (also occurring in 
G. duebeni). The electromorph distribution at 
other polymorphic loci shows similar discontinui¬ 
ties (Table II). 

G. d . duebeni and G. d . celticus are characterized 
by the predominating Apk°, Est-2 a , Got-2 c , Hk f, 
Me c > d , Mptf, Pep-4 C > $, 6Pgd c and Pk c > e electro¬ 
morphs. Substantial differences in electromorph 


distribution between G. d . duebeni and G. d. celti¬ 
cus ^ccur at the Got-1 , Mdh-1 , Me, Pep-4 and Pk 
loci. 

A half matrix of diagnostic loci (0.99 criterion; 
Ayala & Powell, 1972) of area by area comparisons 
of samples of G. gauthieri , G. ibericus and G. due¬ 
beni is summarized in Table III. Values of D, aver¬ 
aged by species and by area, are shown in Table IV. 
Strongly “deviating” samples did not participate in 
this averaging. 

The UPGMA dendrogram (Fig. 2) and the Jelly 
network (Fig. 3) show that levels of intra-specific 
intra-continental genetic differentiation among 
samples of G. gauthieri from within the Iberian 
peninsula (in particular G. gauthieri samples from 
area G with regard to samples from areas E-F) 
may be higher than among inter-continental com¬ 
parisons of samples from North Africa and the 
Iberian peninsula. D-\ alues and the number of di¬ 
agnostic loci of inter-continental comparisons 
among G. gauthieri samples from areas D (stations 
6-9) and G (stations c, d, e) are consistently lower 
than those among G. gauthieri samples from areas 
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+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ 

1.40 1.26 1.12 0.98 0.84 0.70 0.56 0.42 0.28 0.14 0.00 

D 


G. d. duebenf 2 
G. d. duebeni 2 
G. d. ceiticus 3 


G. d. ceiticus 

4 

G. gauthieri 

5 

G. gauthieri 

6 

G. gauthieri 

7 

G. gauthieri 

9 

G. gauthieri 

8 

G. gauthieri 

10 

G. gauthieri 

11 

G. gauthieri 

a 

G. gauthieri 

b 

G. ibericus 

f 

G. gauthieri 

d 

G. gauthieri 

c 

G. gauthieri 

e 

G. ibericus 

9 


Fig . 2. UPGMA dendrogram 
of Nei’s (1972) genetic dis¬ 
tance CD) based upon 20 en¬ 
zyme loci. 


HTU 3 


HTU 5 

HTU 4 N *, 


HTU 6 


HTU 2 


HTU 1 


-C 


-G. gauthieri 7 

G. gauthieri 9 

G. gauthieri 8 


.— G. gauthieri e 

I |- G. gauthieri c 

G. gauthieri d 


G. gauthieri a 
— G. gauthieri b 


■C, 


G. gauthieri 10 
G. gauthieri 11 


G. ibericus f 


G. gauthieri 6 
- G. gauthieri 5 


G. ibericus g 


_I- G.d. ceiticus 4 

_ *~G. d. ceiticus 3 

_|— G. d. duebeni 1 

*— G. d. duebeni 2 

+ — + — +— + — + — + —+— + — + — + —+— + ---+—+— + — + — + — +— +— + — +— +—+— +— + 

0.00 0.16 0.32 0.48 0.64 0.80 0.96 

D 

Fig . 3. Dendrogram of Jelly “character state” network (Ellis, 1987) of 20 enzyme loci, based upon Rogers’ (1984) HAP algorithm, and 
rooted by the outgroup method (Farris, 1972). 






Fig . 4. The Jelly network of Fig. 3, represented as a “Hennigian tree”, with internodia defined by synapomorphous electromorphs. 
* = not strictly contributing to clade definition (plesiomorph undetermined; Patton & Avise, 1983). 


E (stations 10-11) vs. G, and F (stations a, b) vs. 
G, respectively (Tables III and IV). However, inter¬ 
continental ZJ-values among samples from areas D 
and E, and D and F (Table IV) are higher than 
intra-continental values among areas E and F. 
These values are also higher than intra-areal D- 
values of area D, and D-values among station 5 and 
stations 6-9 (Table IV). Inter-continental £>-values 
among G. gauthieri samples from the geographical¬ 
ly closest areas D and E, and the much more distant 
areas D and F, do not differ substantially (Table 
IV). It is also interesting to note that £>-values of 
inter-specific inter-continental comparisons among 
G. gauthieri from area D (stations 7-9) and G. 
ibericus from station f are lower than those among 
G. gauthieri from area E (stations 10-11) and G. 
ibericus from station f (Table IV). 

According to the UPGMA dendrogram (Fig. 2), 
G. gauthieri from North Africa (stations 5-9) 


forms a monophyletic group with the central Iberi¬ 
an G. gauthieri from stations 10-11 (area E), a-b 
(area F) and G. ibericus from station f (area F), thus 
apart from the western Iberian G. gauthieri from 
stations c-e (area G) and G. ibericus from station 
g (area F). Within this cluster, North African popu¬ 
lation samples (stations 5-9) form a clearly distinct 
subcluster, whereas southern Iberian G. gauthieri 
(stations 10-11) clusters with central Iberian G. 
gauthieri + G. ibericus from station f. This is in 
contrast with the Jelly network (Fig. 3), where G. 
gauthieri samples from the North African stations 
7-9, and 5, 6 are represented by successive 
branches between the clusters formed by G. 
gauthieri from area G (western Spain; Fig. 1 A), and 
G. gauthieri + G. ibericus from southern and cen¬ 
tral Spain. Moreover, G. ibericus from stations f 
and g cluster together. 

When the branches of the Jelly network are 
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Table V. Intercontinental values of the mean genetic distance 
( D ) and standard error (s.e.) among stations and areas in North 
Africa and the Iberian peninsula, jack-knifed over 20 enzyme 
loci. Corresponding values of one unit D are listed along; these 
values are based on a separation of 5.5 My; for station legend 
see Table I. 


North Africa 
station(s) area 

Iberian peninsula 
vs. station(s) area 

D 

± 

s.e. 

one unit D 
(My) 

range 

5 


10-11 

E 

0.63 

± 

0.01 

8.73 

0.59-0.70 

5 

- 

a-b 

F 

0.59 

± 

0.01 

9.32 

0.53-0.63 

6 

D 

10-11 

E 

0.54 

± 

0.01 

10.19 

0.48-0.57 

6 

D 

a-b 

F 

0.44 

± 

0.01 

12.50 

0.39-0.47 

7-9 

D 

10-11 

E 

0.41 

± 

0.01 

13.41 

0.35-0.44 

7-9 

D 

a-b 

F 

0.37 

± 

0.01 

14.86 

0.31-0.40 

average 




0.50 



11.50 



defined by synapomorphic electromorphs (Fig. 4), 
internodia defining G. gauthieri from stations c-e, 
and G. ibericus from stations f-g are represented 
by trichotomies. Furthermore, the successive 
branches of the North African stations 7-9, and 5, 
6 form a monophyletic group. The topology of this 
“Hennigian tree” is in outline (except for the 
trichotomies) similar to the UPGMA dendrogram. 

To provide a calibration of the molecular clock, 
inter-continental intra-specific D -\alues among G. 
gauthieri from stations 5, 6, 7-9 (area D), stations 
10, 11 (area E) and a, b (area F) were jack-knifed 
over loci (Table V). According to Hsu et al. (1977), 
the opening of the Strait of Gibraltar occurred 6-5 
My ago. Thus, based on 5.5 My, a Z> mean -value of 
0.50 for inter-continental comparisons corresponds 
to a one unit £>-value of 11.5 My. 


Genetic structure of populations and estimated 
levels of gene flow 

The results of the F-statistic analysis (0) for G. 
gauthieri , G. ibericus , G. d, duebeni and G. d. celti- 
cus are shown in Table VI. Statistically significant 
population subdivision was found for G. gauthieri , 
G. ibericus and G. d . celticus at most polymorphic 
loci: all 16 in G. gauthieri ; all 8 in the two samples 
of G. ibericus , and all 4 in G. d. celticus. This 
degree of subdivision is about the same as recorded 
in combined datasets with G. gauthieri and G. 


Table VI. 0-indices according to Weir & Cockerham (1984) of 
the samples examined. G. gauth . = G. gauthieri; G. ib. = G. 
ibericus; G. d. dueb . = G. d. duebeni; G. d. celt . = G. d. celti¬ 
cus. Significance of population subdivision was tested with a 
chi-square heterogeneity test. * = P < 0.05. Data of Scheep¬ 
maker et al. (1988) are included. 


species 

Area(s) 

Stations 

G. gauth. 
G. ib. 
D-G 

5-11, a-g 

G. gauth. 

D-G 

5-11,a-e 

G. ib. 

F 

f.g 

G. d. dueb. 
G. d. celt. 

-,c 

1-4 

G. d. dueb. 

1,2 

G. d. celt. 

C 

3.4 

locus 

Ada 

1.00* 


1.00* 




Apk 

0.49* 

0.48* 

- 

- 

- 

_ 

Alp-1 

1.00* 

1.00* 

- 

- 

_ 

_ 

Alp-2 

0.89* 

0.87* 

- 

- 

. 

. 

Est-1 

1.00* 

1.00* 

- 

- 

- 

. 

Esl-2 

0.04 

0.88* 

1.00* 

. 

. 

. 

Gdh 

1.00* 

1.00* 

- 

- 

- 

. 

Got-1 

0.99* 

0.86* 

1.00* 

0.65* 

0.05 

0.40* 

Got-2 

0.78* 

0.98* 

0.46* 

0.00 

0.00 

- 

Gpi 

0.65* 

0.66* 

0.59* 

0.59* 

0.05 

0.21* 

Hk 

1.00* 

. 

1.00* 

- 

- 

. 

Mdh-1 

0.91* 

0.90* 

- 

0.87* 

0.02 

0.10* 

Mdh-2 

0.91* 

0.91* 

- 

- 


. 

Me 

1.00* 

1.00* 

- 

0.99* 

0.00 

- 

Mpi 

0.81* 

0.82* 

0.71* 

- 

_ 

- 

Pep-2 

0.91* 

0.88* 

0.99* 

- 

. 

. 

Pep-4 

0.86* 

0.80* 

- 

0.32* 

0.04* 

0.36* 

Pk 

0.94* 

0.94* 

- 

1 . 00 * 

- 

- 

Mean 

0.89 

0.89 

0.89 

0.74 

0.04 

0.29 


ibericus (Tables VI, VII). However, much less 
population subdivision was found in the brack¬ 
ish/marine form G. d . duebeni . Creating subsets of 
samples of G. gauthieri by omitting the most 
“deviating” samples, by combining samples by 
area, or by pairwise combinations of geographical¬ 
ly close samples does not change this pattern. These 
subsets show that statistically significant popula¬ 
tion subdivision still exists over small geographic 
ranges (Table VII). 

Mean values of 0 within species and subsets of 
population samples are listed in Tables VIII and 
IX, along with the derived Am-values (Nm e ). The 
average frequencies of private alleles p (1) and the 
derived Nm - and Am*-values (Slatkin, 1985) were 
calculated for sets of four samples and more. 

Genetic drift and bottlenecks will result in sub¬ 
stantial local differentiation if Nm < 1; when Nm 
> 1, populations tend to become panmictic for 
neutral alleles (e.g. Wright, 1931). Except for G. d. 
duebeni (stations 1 and 2) and G. gauthieri from 
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0.88 0.87 0.87 0.84 0.80 0.73 0.77 0.69 0.72 0.83 0.70 0.09 0.38 0.50 0.39 0.26 
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Table VIII. Various Nm estimates in G. gauthieri , G. ibericus , G. d. duebeni and G. d. celticus. n.p.a. = number of private alleles; 
p (1) = average frequency of private alleles (Slatkin, 1985); Nm = estimate derived from p (1); Nm * = idem, with correction for sample 
size; 6 = /^statistic (Weir & Cockerham, 1984); Nm d = estimate derived from 0; km (av) = average geographic distance among sam¬ 
ples in kilometres. Data of Scheepmaker et al. (1988) are included. 


Species 

Area(s) 

Stations 

No. of 
samples 

0 

n.p.a. 

m 

Nav 

Nm 

Nm* 

Nm 0 

km(av) 

G. gauthieri 

D,E,F,G 

5-11; a-g 

14 

0.89 

15 

0.58 

41.08 

0.02 

0.01 

0.03 

403 

+ G. ibericus 

G. gauthieri 

D,E,F,G 

5-11; a-e 

12 

0.88 

13 

0.47 

39.20 

0.04 

0.02 

0.04 

423 

G. gauthieri 

D,E,F 

5-11; a-b, f-g 

11 

0.89 

19 

0.50 

38.30 

0.03 

0.02 

0.03 

380 

+ G. ibericus 

G. gauthieri 

D,E,F 

5-11; a-b 

9 

0.87 

15 

0.42 

37.58 

0.04 

0.03 

0.04 

405 

G. gauthieri 

D,E 

5-11 

7 

0.87 

14 

0.43 

35.71 

0.04 

0.03 

0.04 

256 

G. gauthieri 

D 

5-9 

5 

0.84 

15 

0.42 

35.80 

0.04 

0.03 

0.05 

235 

G. gauthieri 

D 

6-9 

4 

0.80 

15 

0.58 

35.82 

0.02 

0.02 

0.06 

34 

G. gauthieri 

D 

7-9 

3 

0.73 

- 

- 

- 

- 

- 

0.09 

41 

G. gauthieri 

E, F 

10-11; a-b 

4 

0.83 

7 

0.55 

39.80 

0.03 

0.02 

0.05 

400 

G. gauthieri 

G 

c-e 

3 

0.38 

- 

- 

- 

- 

- 

0.41 

40 

G. d. duebeni 

-,C 

1-4 

4 

0.74 

4 

0.05 

36.03 

3.13 

2.17 

0.09 

341 

+ G. d. celticus 


area E (stations 10 and 11), in all species and subsets 
examined, Am-values derived from 0 ( Nm e ) are 
less than one (Tables VII, IX), suggesting an impor¬ 
tant degree of genetic isolation. Except for the G. 
gauthieri samples from stations 10 and 11, there is 
no correlation between Aw^-values and geographic 
distance (Tables VIII and IX). These data suggest 
that substantial barriers to gene flow occur over 
small geographic ranges. 

Ata-estimates according to Slatkin, 1985 (in par¬ 
ticular estimates uncorrected for sample size) show 
a very close agreement with Afo^-estimates (Table 
VIII); however, the two estimates contrast in the 
dataset formed by G. d . duebeni and G. d. celticus 
(stations 1-4). 

The results of the CAF analysis (Slatkin, 1981) of 
gene flow levels in different sets and subsets of G. 
gauthieri populations (whether or not in combina¬ 
tion with the related form G. ibericus) are repre¬ 
sented in Figs. 5 and 6. All gene flow profiles show 
high initial p (/)-values, and almost all exhibit wide 
fluctuations. These curves resemble closely those 
presented by Slatkin (1981) as characteristic for 
low-gene-flow species (see Scheepmaker, 1990). 


Table IX. Various Nm estimates in pairwise combinations of 
samples of G. gauthieri , G. ibericus , G. d. duebeni and G. d. cel¬ 
ticus. For legend see Table VIII. 


Species 

Area 

Stations 

9 

Nm e 

km(av) 

G. gauthieri 

G 

c, e 

0.39 

0.39 

48 

G. gauthieri 

G 

d, e 

0.50 

0.25 

54 

G. gauthieri 

G 

c, d 

0.26 

0.72 

17 

G. gauthieri 

E 

10,11 

0.09 

2.38 

8 

G. gauthieri 

D 

7,9 

0.77 

0.08 

44 

G. gauthieri 

D 

8,9 

0.69 

0.11 

19 

G. gauthieri 

D 

7,8 

0.72 

0.10 

64 

G, gauthieri 

F 

a, b 

0.70 

0.11 

35 

G. ibericus 

F 

f, g 

0.89 

0.03 

26 

G. d. duebeni 

- 

1,2 

0.04 

5.98 

400 

G. d. celticus 

C 

3,4 

0.29 

0.61 

24 


Cross-breeding experiments 

The results of the cross-breeding experiments are 
listed in Table X. Due to high mortality rates (prob¬ 
ably caused by unknown compounds in the Amster¬ 
dam tap water) in population samples from 
southern Spain (area E) and Morocco (area D), ex¬ 
perimental conditions were not always optimal. All 
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Table X. Results of cross-breeding experiments between samples of G. gauthieri and G. ibericus from North Africa and the Iberian 
peninsula. g 6 = G. gauthieri from station 6; g J0 = G. gauthieri from station 10; g b = G. gauthieri from station b; i g = G. ibericus 
from station g; for legend stations see Table I. 


Cross 

99 x crcr 

N of 9 9 

N of crcr 

N of 

ovigerous 9 9 

N of ovigerous 

9 9 with offspring 

F2 

obtained 

s 6 

x 8 6 

122 

84 

53 

32 

+ 

810 

x 8io 

7 

6 

3 

1 

- 

8b 

X 8b 

17 

17 

14 

8 

+ 

ig 

* 'e 

50 

50 

30 

10 

+ 

8s 

x iio 

60 

16 

- 

- 

- 

8 io 

X 86 

23 

26 

8 

1 

- 

86 

X ig 

36 

36 

13 

9 

+ 

ig 

X 8 6 

84 

52 

4 

- 

- 

86 

X 8 b 

50 

31 

18 

4 

- 

8b 

X 8 6 

107 

90 

6 

- 

— 


Table XI. Distribution of numbers of electromorphs common in western European, Iberian and North African members of the G. pulex- 
group in samples from North Africa and the Iberian peninsula; for legend stations see Table I. 


species 

G. gauthieri 

G. gauthieri 

G. gauthieri 

G. gauthieri 

G. gauthieri 

G. ibericus 

G. ibericus 

G. gauthieri 

area 

- 

D 

D 

E 

F 

F 

F 

G 

station (s) 

5 

6 

7-9 

10-11 

a, b 

f 

g 

c-e 


Morocco 

Morocco 

Morocco 

Spain 

Spain 

Spain 

Spain 

Spain 

Common in western European 
members of the G. pulex- group 

15 

18 

27 

15 

16 

20 

17 

17 

Common in G. duebeni 

8 

6 

11 

7 

8 

8 

11 

10 

Common in western European members 
of the G. pulex- group; 
present in Iberian samples; 
lacking in North African samples 




2 

1 

4 

6 

4 

Common in western European members 
of the G. pulex- group; 
present in North African samples; 
lacking in Iberian samples 

2 

2 

5 






Restricted to both North African and 
Iberian samples 

5 

5 

8 

4 

3 

3 

3 

6 

Exclusively in North African samples 

0 

1 

0 

- 

- 

- 

- 

- 

Exclusively in Iberian samples 

- 

- 

- 

2 

2 

0 

3 

2 


reference crosses but one (g I0 X g }0 ) yielded F2 
offspring. In three combinations, F2 could not be 
obtained in reciprocal crosses: g 6 x g I0 , i g x g 6 , 
and g b x g 6 . In one cross, g 6 x i g9 F2 could be ob¬ 
tained (g refers to G. gauthieri ; i to G. ibericus ; 6, 
10, b, and g to stations 6,10, b, and g, respectively). 


Discussion and conclusions 

Electromorph distribution with reference to central 
Iberian G. gauthieri samples and related taxa of 
western Europe north of the Pyrenees 

Differences in electromorph distribution among 
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G. gauthieri, G. ibericus 

(Stations 5-11; a-g ) 

— 

— 

G. gauthieri 

(Stations 5-11) 

G. gauthieri 

(Stations 5-11; a-e) 

— 


G. gauthieri 

(Stations 5-9) 

G. gauthieri, G. ibericus 

(Stations 5-11; a-b, f*g) 

◄— 

—► 

G. gauthieri 

(Stations 5-8) 

G. gauthieri 

(Stations 5-11; a-b) 


— 

G. gauthieri 

(Stations 10-11. a-b) 




Fig. 5. CAF gene flow profiles (Slatkin, 1981) in G. gauthieri 
and G. ibericus. p (i) = conditional average frequency (CAF) of 
an electromorph; i/d — incidence of an electromorph, where 
/ = the total number of samples it appears in, and d = the total 
number of samples examined. The points plotted are average 
values of p (/) for electromorphs having the same i/d-v alue. For 
station legend see Table I. 

southern Iberian and North African G. gauthieri 
samples are less conspicuous than among samples 
within the Iberian peninsula; however, with regard 
to taxa north of the Pyrenees, G . gauthieri is much 
more differentiated. Electromorph frequency data 
of G. gauthieri , G. ibericus and members of the G. 
pulex- group in northwestern Europe were reported 
by Scheepmaker et al. (1988), Scheepmaker & Van 
Dalfsen (1989) and Scheepmaker (1990), respective¬ 
ly. The numbers of electromorphs in population 
samples of G. gauthieri and G. ibericus from North 
Africa and the Iberian peninsula in common with 
members of the G. pw/ex-group in North-West 
Europe are listed in Table XI. The highest number 
of “western European’’ electromorphs was found 
in area D (North Africa); however, samples from 


Fig. 6. CAF gene flow profiles in subsets of G. gauthieri . For 
legend see Fig. 5. 

this area also exhibit the highest levels of polymor¬ 
phism. The combination of the fixed electromorphs 
Alp-l b and Alp-2 C for instance, that occurs in all 
North African samples, is typical for all western 
European taxa of the G. pulex- group (this combi¬ 
nation, also found in the outgroup G. duebeni, is to 
be considered plesiomorphic). The Gpi e 'f>*> h > 1 
and the Mpi h > 1 electromorphs are common in 
western European members of the G. pulex- group 
and North African G. gauthieri , but lack in most 
populations of G. gauthieri from the Iberian penin¬ 
sula, etc. The simultaneous occurrence of fixed 
Apk° t 6Pgd h and Pk° electromorphs (occurring in 
all North African samples) is never found in mem¬ 
bers of the G. pulex- group in North-West Europe 
and characterizes G. gauthieri (and G. ibericus) 
from the Iberian peninsula; likewise, the Est-2 b 
electromorph is fixed in all Iberian and North Afri¬ 
can samples of G. gauthieri and is rare among 
western European members of the G. pulex- 
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group. Although North African G. gauthieri popu¬ 
lations exhibit some (in part plesiomorphic) charac¬ 
ters that are typical for North-West European 
members of the G. pulex- group, electromorphs 
characterizing taxa from the Iberian peninsula 
dominate. 

The Apk a , Est-2 b and 6Pgd b electromorphs are 
shared by all samples of G. gauthieri from area G 
(stations c-e), and G. gauthieri and G. ibericus 
from areas D-F (stations 5—11, a-b and f-g). 
However, as the plesiomorph remains undeter¬ 
mined, they do not contribute to clade definition 
(although they may provide supporting evidence 
for previously defined clades; Patton & Avise, 
1983). Therefore, no clade can be defined including 
all these samples (Fig. 4). 

Electromorphs occurring only in either all Iberi¬ 
an or all Morroccan population samples of G. 
gauthieri were not found. On the other hand, the 
highest number of ‘‘endemic” electromorphs were 
found among G. gauthieri plus G. ibericus from the 
Iberian peninsula. Moreover, inter-continental 
inter-areal D-v alues among areas D and G are lower 
than intra-continental inter-areal £>-values among 
areas E and F vs. area G. These data suggest that 
not all populations accumulate evolutionary 
change at the same rate (i.e., faster rates occur in 
taxa from the Iberian areas E and F). 


Genetic differentiation , origin and dispersal 

According to Fig. 2, western Iberian G. gauthieri 
from area G cluster apart from North African, and 
southern and central Iberian G. gauthieri from sta¬ 
tion 5 and areas D, E and F. North African G. 
gauthieri samples (stations 5-9) form a clearly dis¬ 
tinct subcluster among southern and central Iberian 
G. gauthieri and G. ibericus samples from areas E 
and F. This pattern suggests that North African 
populations of G. gauthieri originate from the 
Iberian peninsula, in particular from common an¬ 
cestors with G. gauthieri from areas E and F. 

According to the Jelly network (Fig. 3), a western 
European ancestor gave rise to (ancestors of) G. 
duebeni , and next, to common ancestors of Iberian 
and North African G. gauthieri (HTU 1 in Fig. 3). 


North African G. gauthieri populations (station 5, 
area D) originate in 4 successive branches from 
common ancestors with G. gauthieri (and G. iberi¬ 
cus) from areas E and F (HTU’s 2-6). The 4 succes¬ 
sive branches of stations 5-6, 7, 8 and 9 suggest 
that the Strait of Gibraltar was traversed four times 
(or a commuting movement). However, if these 
branches (HTU’s 2-6) are considered as a polyto- 
my, or as a monophyletic group (as in the “Hen- 
nigian tree” of Fig. 4), traversing once suffices. 

From HTU 1 to HTU 2, important changes in 
frequency distribution occur at the Gpi , Pep-2 and 
Pep-4 loci. However, HTU 2-6 (i.e. the 4 succes¬ 
sive branches) are rather similar and differ essen¬ 
tially by gradual changes at the Gpi , Got-l y Me and 
Mpi loci (except for HTU 3 to HTU 4, where the 
plesiomorphic fixed Gdh b electromorph is replaced 
by Gdh a ). Therefore, it may be justified to repre¬ 
sent HTU 2-6 as a polytomy. As for the remaining 
loci, HTU 1-6 are all characterized by the fixed 
Alp-l b y Alp-2 c y Est-l b (all plesiomorphic), Est-2 b y 
Hk b y and 6Pgd a y and predominant Apk° and Pk° 
electromorphs (for sake of space, HTU frequency 
data are not reproduced here). 

As pointed out by Scheepmaker et al. (1988), 
Iberian and North African G. gauthieri and G. 
ibericus are likely to originate from northwestern 
European and not from African ancestors. Gam- 
marids are lacking in the rest of Africa, and a major 
adaptive shift of a species group from temperate to 
more extreme conditions seems improbable (Goli¬ 
kov & Tzvetkova, 1972). Unless the extension of 
(ancestors of) G. gauthieri included the Tyrrhenian 
and Calabrian areas (which are at present devoid of 
G. gauthieri- like forms; Pinkster, pers. comm.), 
and unless G. gauthieri reached North Africa from 
the East after the junction of the Eurasian and Afri¬ 
can plates, the most obvious dispersal route is via 
the northeastern part of the Iberian peninsula (be¬ 
fore the upheaval of the Pyrenees), and from there 
to southern Spain and North Africa. 

Geographical barriers to gene flow and vicariant 
events 

Today, geographic and climatological factors pre¬ 
vent the genetic exchange among areas D and G. 
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Area G is at present primarily separated from the 
other areas by the Hercynian Sierra de Gredos, 
Sierra de Guadarrama and Celtiberic Cordillera; 
this barrier may have been effective since pre- 
Messinian times (Ager, 1980; Scheepmaker et al., 
1988). Areas F and E are separated by the more re¬ 
cent, Alpine Betic Cordillera, and areas D and E-F 
by the Strait of Gibraltar. 

In addition to recent and historical barriers, pre- 
Messinian subdivision of ancestral populations 
may also account for the genetic divergence record¬ 
ed among G. gauthieri and G. ibericus from areas 
D, E and F. Since the Eocene-Oligocene boundary, 
the emerged eastern part of the Iberian peninsula 
(reaching from the Betic Cordillera via the Balearic 
promontory, Sardinia and Corsica up to southern 
France) was separated from the western part by ma¬ 
rine transgressions (e.g. the formation of the Valen¬ 
cia trough) and fragmented repeatedly by tectonic 
activity and successive marine transgressions and 
retreats (Pomerol, 1973; Dercourtet al., 1985; Mal¬ 
donado, 1985). The Betics (including area E) must 
have resembled an archipelago during several 
phases of its history (Notenboom, 1988). 

Ancestors of G. gauthieri were possibly already 
present on the former Betic Archipelago (including 
area E and the Moroccan Rif area) and the formerly 
emerged eastern part of the Iberian peninsula, as 
probably were (ancestors of) Echinogammarus 
klaptoczi ebusitanus Schaferna, 1908 and E. sicilia- 
nus sicilianus Karaman & Tibaldi, 1973. The latter 
two species are actually known from Mallorca, 
Menorca, Ibiza, Malta and Libya (North Africa), 
and from Mallorca and Sicily, respectively (Stock, 
1978; Stock, pers. comm.). However, the absence 
of G. gauthieri or any related forms in other rudi¬ 
ments of this former land such as the Balearic 
Islands and Sardinia (Pinkster, pers. comm.) sug¬ 
gests that (an ancestor of) G. gauthieri has not ex¬ 
tended to this area until the Messinian salinity crisis 
(6-5 My ago). Unless G. gauthieri invaded North 
Africa from the East after the junction of the Eura¬ 
sian and African plates, this was likely to happen 
from the Iberian peninsula. 

/7-values among samples from areas E and F are 
not much higher than among North African and 
Iberian samples from areas E and F; /7-values 


among North African and Iberian samples from 
area G are even lower than among samples from 
areas E-F, and area G. If G. gauthieri did reach 
North Africa from the East, these values would 
have been (presumably) much higher. 

In North Africa, /7-values of samples from sta¬ 
tion 6 vs. stations 7-9 (all area D; mean distance c. 
50 km) are higher than those between stations 5 and 
6 (High Atlas and Middle Atlas respectively; geo¬ 
graphic distance c. 420 km; Tables IV, VIII). As 
these areas were separated by marine transgressions 
during the Pliocene, the presence of G. gauthieri in 
the High Atlas may be relatively recent. 

Today, desertification as a result of deforestation 
and erosion in many areas (e.g. central Spain, parts 
of the Betic Cordillera, the Rif zone in North Afri¬ 
ca, etc.) represents important barriers among some 
of these areas. Due to semi-arid conditions and pol¬ 
lution gammarids disappeared from many locali¬ 
ties. The resulting degree of discontinuity among 
conspecific samples (even within the same area or 
drainage system) makes it difficult to find relation¬ 
ships among areas. However, as intra-continental 
/7-values among both North African and Iberian 
population samples are lower than inter-conti¬ 
nental values among North Africa and the Iberian 
peninsula, the opening of the Strait of Gibraltar is 
probably the most important among factors limit¬ 
ing gene flow. 


Calibration of the molecular clock 

An initial calibration of the molecular clock can be 
performed when inter-continental values of D ex¬ 
ceed intra-continental /7-values. Samples of G. 
gauthieri from areas E and F (forming a monophy- 
letic group in Figs. 2 and 3) satisfy this requirement. 
However, the data set reported here represents a 
serious risk of pitfalls. The range of /7-values 
resulting from the jack-knife procedure is 
0.31-0.70 (Table V). Moreover, the opening of the 
Strait of Gibraltar may be dated at 5, 5.5 or 6 My 
ago. Consequently, one unit D may range from 7.1 
to 19.3 My based on the same data set. In addition, 
the stations sampled do not necessarily represent 
the geographically closest areas on either side of the 
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Strait of Gibraltar; genetically closer forms may 
have occurred in areas being devoid of gammarids 
today (e.g. the Rif zone and the southernmost 
Betics). 

Rates of accumulation of genetic change may be 
strongly biased by stochastic processes. Averaging 
over several samples from an area may not compen¬ 
sate for this effect. Thus, the present calibration 
should not be applied to species that are much less 
subjected to stochastic processes (e.g. G. d . duebe- 
ni , G. zaddachi and G. tigrinus; Scheepmaker, 
1990). Moreover, populations on either side of the 
Strait of Gibraltar possibly do not accumulate 
evolutionary change at the same rate (e.g. areas E, 
F). Consequent to these uncertainties, other data 
sets and well-dated geological events are required to 
provide conclusive evidence for the meaning of the 
present calibration. With reference to literature 
data, a one unit D -\alue of 11.5 fits in the range of 
values (5-20 My) proposed by different authors 
(e.g. Nei, 1975; Sbordoni et al., 1980; Thorpe, 
1982; Bert, 1986). 


Genetic differentiation and estimated levels of gene 
flow 

Populations of G. gauthieri and G. ibericus are 
highly subdivided at many loci. Single locus 6- 
values of one are no exception, and one may 
wonder if the samples examined are still to be con¬ 
sidered conspecific. Significant population subdivi¬ 
sion occurs at much higher levels than those record¬ 
ed in cave-dwelling species (Caccone & Sbordoni, 
1987) or in related members of the G. pulex-g roup 
from central and southern France (Scheepmaker, 
1990). Estimated levels of gene flow (Nm*, Nm and 
Nm e ) in sets, subsets, or pairwise combinations of 
populations are in all but one case (G. gauthieri 
from stations 10-11) much smaller than one. This 
is in contrast with results of Scheepmaker (1990), 
who found occasionally high Nm*, Nm and Nm e - 
levels in population subsets and pairwise combina¬ 
tions. As pointed out by Scheepmaker (1987,1990), 
periodic desiccation of riverbeds and erosion by 
heavy rainfall may be the acting force of population 
fragmentation and consequent stochastic processes 
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in the Mediterranean area. The degree of tem¬ 
porary isolation of population fragments subjected 
to genetic drift increases with the xerothermic index 
(Bagnouls & Gaussen, 1953; Scheepmaker, 1987) 
and the length of the dry season. 

The xerothermic index of the stations sampled 
for this study ranges from about 70 (Malaga) to 93 
(Amizmiz, 7 km from station 5 and 30 km from 
Marrakech). However, sampling areas are often 
separated by areas with a much higher xerothermic 
index; in Marrakech for instance, this value is 178, 
and 118 in Fes (Fig. 1A, near area D). The mean 
length of the dry season is 7 months in Marrakech 
(Bagnouls & Gaussen, 1953); in periods of drought 
this may be much longer. The consequent semi¬ 
permanent desiccation of riverbeds and water¬ 
courses may result in long-term isolation of popula¬ 
tions and increased levels of genetic divergence. 
During the Pliocene-Pleistocene sequence of cold 
phases similar dry conditions resulting in popula¬ 
tion subdivision may have prevailed (Scheepmaker 
et al., 1988), whereas intermittently more temperate 
pluvious conditions during the Pleistocene intergla¬ 
cials (Biberson, 1970) possibly allowed range ex¬ 
tension. 


Cross-breeding experiments 

In cross-breeding experiments, inter-continental G. 
gauthieri populations from stations 6, 10 and b 
(D = 0.41 -0.53) suggest partial intersterility; F2 
was never obtained. However, negative results may 
be due to unfavourable laboratory conditions; 
Scheepmaker et al. (1988) obtained F2 juveniles in 
one of the reciprocal crosses between genetically 
more distant G. gauthieri samples (Z> = 0.60). 
Moreover, in inter-specific inter-continental cross¬ 
es between G. gauthieri and G. ibericus from sta¬ 
tions 6 and g (D = 0.82), one of the reciprocal cross¬ 
es yielded F2 offspring. These results agree with 
findings of Scheepmaker et al. (1988), and cor¬ 
roborate their opinion that cross-breeding experi¬ 
ments are not conclusive in attempts to resolve the 
taxonomic status of Iberian members of the G. 
pulex- group. 

The results of the cross-breeding experiments 
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(i.e. at least partial interfertility among genetically 
distant population samples) suggest that, although 
important differences occur in structural genes cod¬ 
ing for enzymes, regulatory genes may still be com¬ 
patible to a high degree. Most regulatory genes are 
active during the ontogenetic process. To complete 
this process, compatibility of regulatory genes is 
required (Dubois, 1983). The conservatism of regu¬ 
latory genes is supported by the morphological 
homogeneity among genetically distant popula¬ 
tions of G. gauthieri (and, to some extent, G. 
ibericus). 
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